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Abstract 

We present a numerical study of the effect of the repulsive logarithmic inter-particle interaction on 
the ground state configuration and the frequency spectrum of a confined classical two-dimensional 
cluster containing a finite number of particles. In the case of a hard wall confinement all particles 
form one ring situated at the boundary of the potential. For a general r" confinement potential, 
also inner rings can form and we find that all frequencies lie below the frequency of a particular 
mode, namely the breathing-like mode. An interesting situation arises for the parabolic confined 
system (i.e. n = 2). In this case the frequency of the breathing mode is independent of the 
number of particles leading to an upper bound for all frequencies. All results can be understood 
from Earnshaw's theorem in two dimensions. In order to check the sensitivity of these results, the 
spectrum of vortices in a type II superconductor which, in the limit of large penetration depths, 
interact through a logarithmic potential, is investigated. 
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I. INTRODUCTION 



In recent years, classical clusters with a finite number of particles moving in a two- 
dimensional plane were studied for different kinds of confinement and interaction potentials. 
Bedanov et al. studied the classical system with a finite number of particles interacting 
through a repulsive 1/r potential and moving in a two-dimensional plane, confined by a 
parabolic potential. The 1/r potential is the Coulomb potential in a three-dimensional 
world. This theoretical system models experimental realizations, such as electrons on the 
surface of liquid helium ^l? electrons in quantum dots in high magnetic fields ^, colloidal 
suspensions 5| and confined plasma crystals j^. A detailed investigation of the structure l| 
and the spectral properties 2J of these 2D clusters was carried out. Recently, ground state 
and metastable configurations and/or spectral properties of 2D clusters with a parabolic 
confinement potential but different interaction potentials were investigated (see Ref. 
0, y, y, 13 for a I/'" potential , R ef. lol 11 1 for a dipole interaction, Ref. [l^ [ll| for a 
logarithmic potential, Ref. [lO, ll^ for screened Coulomb interaction, ...) and with a 1/r 
Coulomb interaction potential but with different confinement potentials (see Ref. for a 
Coulomb confinement potential, Ref. Q| for a hard wall potential). 

Laughhn has shown that for strongly interacting electrons in a magnetic field at 
certain fractional fillings (1/3, 1/5, . . .) of the Landau levels, the maxima in the pair correla- 
tion function correspond to the equilibrium positions in a classical one dimensional plasma 
where the electrons interact via a repulsive logarithmic potential. The idea was extended 

n y ^ y l_. 

to quantum dots [16|| with few electrons as well as bosons and to rings jl7( with a few elec- 
trons. It was shown that the electrons actually get localized at the positions where the pair 
correlation function peaks and the excitations of the system can be understood as rotational 
modes of the center of mass and vibrational modes. The eigenenergies of the system can 
be easily determined from the frequencies of the classical normal modes to a high degree of 
accuracy. 

In this paper we consider systems with various confinement potentials containing particles 
interacting through a repulsive logarithmic interaction. As the solution of the 2D Poisson 
equation yields a logarithmic potential, one can say that the Coulomb potential in a 2D 
world is logarithmic. Vortices in a film of liquid Helium interact through a logarithmic 
potential [l^. Also vortices in a type II superconducting 2D film, for a low concentration 
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of vortices, is expected to interact through a logarithmic potential The logarithmic 

interaction between vortices was used to study the stable vortex configurations in a disk 



shaped superconductor [20|. Since the vortices are infinitely long in the 2;-direction, the 
symmetry makes it effectively 2D. Hence the model may have some experimental realizations 
and it may be possible to verify experimentally the predictions made here. 

In Section II we describe our model systems. The results for the ground state config- 
urations are shown in Section III, and for the spectra in Section IV. Section V tests the 
sensitivity of the results by looking at the spectrum of vortices in a type II superconductor, 
which interact through a logarithmic potential for a large penetration depth. Our conclu- 
sions are given in Section V. 

II. MODEL AND NUMERICAL APPROACH 

The Hamitonian of a 2D system of charged particles in a r" confinement potential and 
interacting through a repulsive logarithmic potential is given by 

1=1 i>j 

where m is the mass of the particle, ujo the radial confinement frequency, and fi = {xi, Ui) the 
position of the ith particle with = |rj|. A hard wall confinement is obtained for n — oo 
and in this case R equals the radius of the hard wall. The interaction potential is taken to 
be logarithmic: 

V{r) = -p\Yi{r/R). (2) 

We can write the Hamiltonian in a dimensionless form if we express the coordinates 
and energy in the following units: r' = P^/"'a~^^"'R^^~'^^^"' , E' = P, with a = ^mujQ. The 
dimensionless Hamiltonian is given by 

TV N 
i=l i>j 

In the limit of a hard wall confinement, the lengthunit becomes r' R. 

We also present results for the actual interaction potential for vortices in a type II su- 
perconductor for A ^ oo as well as for smaller values of A (with A the penetration depth). 
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Vortices in a type II superconductor interact with a potential \21 \ 



VJr) = PKo 



A 



(4) 



with Kq the zero-order Hankel function with imaginary argument, f3 = ^q/(27tX'^) and 
$0 = hc/2e the flux quantum. In the limit of A ^ oo, this potential becomes logarithmic: 



Vsi^) — > — /31nr, for A — > oo (up to a constant). 
With the above units, the dimensionless form of this interaction potential is 

K(r) = Ko 



(5) 



(6) 



where now also A is expressed in the unit r'. 

The numerical method to obtain the ground state configuration is based on the Monte 
Carlo simulation technique supplemented with the Newton method in order to increase the 
accuracy of the energy. The latter technique is outlined and compared with the Monte Carlo 



technique in Ref. 
dynamical matrix 



The eigenmode frequencies are obtained from the eigenvalues of the 



where {r^ j} is the ground state configuration. The eigenvalues of the dynamical matrix are 
the squared eigenfrequencies of the system. The eigenfrequencies are expressed in the unit 
Lu' = ^E'/r'^/m. 



III. GROUND STATE CONFIGURATIONS 



In this section we will discuss the numerically obtained ground state configurations for 
different powers of the confinement potential and show how they can be understood from 
Earnshaw's theorem. Fig. ^ shows as an example the ground state configurations for 40 
particles for (a) a hard wall confinement (i.e. n —>■ oo), (b) a confinement with n = 3 and 
finally (c) a parabolic confinement (i.e. n = 2). One can see that in the case of hard wall 
confinement all particles are situated on one ring at the boundary of the potential. This is 
true for any number of particles in the system, also for clusters with many particles. When 
the confinement differs from the hard wall confinement, particles can be situated in the 
central region. Notice that the density of particles for n = 3 is much larger at the edge than 
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(a) 



(b) 



(c) • 



FIG. 1: Ground state configurations for 40 particles interacting through a logarithmic potential 
for (a) a hard wall confinement, (b) a confinement with n = 3 and (c) for a parabolic confinement. 
The scale is different in each figure, but the distance between the ticks is always one length unit. 
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FIG. 2: The eigenfrequencies for 2 up to 50 particles interacting through a logarithmic potential 
for (a) a n = 3 confinement potential and (b) a parabolic confinement potential. 



;he density is uniform. These 



22 1 in two dimensions. Since in 



in the central region, while for the parabolic confinement 
observations can be understood from Earnshaw's theorem 
the present case the inter-particle potential is logarithmic we can apply the 2D version of 
Gauss's theorem 

/ E • rfr ~ [ p{r)dS, (8) 
Jc Js 

where S is the surface enclosed by C, and p{r) is the enclosed charge. Now it is easy to 



show that charged particles on a 2D surface can not have a stable static equilibrium when 
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FIG. 3: (a) The breathing-like mode for = 40 corresponding to the highest frequency for n = 3. 
The breathing modes for = 40 in the case of a parabohc confinement. 

they interact through a logarithmic potential (i.e. the 2D Coulomb potential). If we take a 
small circle anywhere in the 2D region with no enclosed charges, then, since the line integral 
of the field along the circle has to be zero, the number of outgoing lines of force are equal 
to the number of incoming ones. As long as there are outgoing lines of force, a charge 
left inside this circle can lower its energy by moving along it. Thus in case of hard wall 
confinement, particles can not be in a stable static equilibrium if they are not in contact 
with the hard wall potential. Our simulations indeed show that all particles are pushed on 
one ring situated at the boundary of the hard wall. This is different from the system with 
a hard wall confinement in which the particles interact through a 1/r potential. In Ref. 3] 
it was shown that stable multiple ring structures can be found if a sufficiently large number 
of particles are present in the 2D system. 

As shown above the logarithmic interaction potential alone can not result in a stable 
multi-ring configuration. It is only possible if the hard wall confinement is changed into 
some soft wall r" confinement potential (as shown in Figs, ^b) and (c)). In this case, 
when a particle in the central region is moved in the outward radial direction, the restoring 
force is given by the confinement potential alone. This must balance the nonrestoring forces 
provided by the Coulomb repulsion with the other particles. For the parabolic confinement 
case, this results in a uniform density as the restoring force is everywhere the same. 

IV. THE EIGENMODE SPECTRUM 

Also the eigenmode spectrum for systems with a logarithmic interaction potential differs 
from all the previously studied cases. Again the results can be understood from Earnshaw's 
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theorem. The eigenfrequencies of the system with logarithmic interaction for a n = 3 
confinement potential are, as an example for a soft wall confinement, shown in Fig. |2fa) for 
= 2 up to 50. Notice that there is a pronounced highest frequency branch. We checked 
the corresponding eigenmodes and they all show the same behaviour: all outer particles 
move in the radial direction and have a large amplitude. An example of these breathing-like 
modes is shown in Fig. El^a) for = 40. It is interesting to look now at the eigenfrequencies 
in the case of a parabolic confinement potential, which is shown in Fig. Efb). Again the 
mode corresponding to the highest frequency is the breathing mode (shown in Fig. Efb) for 
A^ = 40). But, as the frequency of the breathing mode in the parabolic case is independent 
of the number of particles (and given by cubreathing = i^max = 2 |23|J), there is an upper bound 
for the frequencies. 

For confinement potentials with n > 2 the potential is steeper than the parabolic con- 
finement at large distances from the center. So for the confinement with n > 2 a mode with 
larger velocities at the boundary than at the center overtakes the breathing mode to be the 
highest mode. But in all these modes there is a simultaneous contraction of a large number 
of particles so we can speak of breathing-like modes. 

Normal modes are effectively a motion of free particles moving in a resultant potential 
created by the other particles and the confinement potential. So whenever we try to move 
a particle from its static equilibrium position, there are restoring forces that bring it back. 
Then, addition of extra particles complicates the effective potential in which the individual 
particles are moving and creates higher normal modes. This potential can become steeper 
and steeper as more particles are added. This results in higher and higher normal modes. 
But this argument breaks down for the logarithmic interaction for which Gauss's theorem 
can be applied: as already shown in previous section, the confinement potential is the only 
restoring force in the outward radial direction and forms an upper bound to the restoring 
force. Addition of new particles does not create a more complicated effective potential but 
it can only result in a slight increase in the frequency of the breathing-like mode due to the 
increase of the nonrestoring Coulomb force on the outer particles. 
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FIG. 4: The spectrum for the paraboUically confiend system interacting through the potential 
Vs{r) as a function of the number of particles for a) A = 1, b) A = 10 and c) A = 50. 




FIG. 5: The eigenmodes corresponding to highest frequency for = 10 for a) A = 1, b) A = 10 
and c) A = 50. 

V. VORTICES IN SUPERCONDUCTORS 



As mentioned in the Introduction, the logarithmic potential is sometimes a good model for 
the interaction between vortices in a type II superconductor. In order to check how important 
small deviations from this logarithmic potential are, we also studied the spectrum of vortices 



in a type II superconductor for a more generally valid interaction potential. In Ref. 
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FIG. 6: The eigenmodes corresponding to highest frequency for = 20 for a) A = 1, b) A = 10 
and c) A = 50, and corresponding to the highest but one frequency for d) A = 1, e) A = 10 and f) 
A = 50. 

is shown that the interaction potential between vortices in a type II superconductor is given 
by Eq. if the penetration length A ^ the coherence length. In the limit of A — oo, this 
potential reduces to the logarithmic potential studied in the previous sections. 

The spectrum for the parabolically confined system is shown in Fig |3] for penetration 
lengths A = 1, 10 and 50. One can clearly see that for A = 1 no upper bound frequency is 
found, while it is for A = 50. That the logarithmic limit is not at all reached for A = 1 is 
also seen in Fig. which shows the eigenmode corresponding to the largest frequency for 

= 10. As shown above, this eigenmode corresponds with the breathing mode in case of a 
logarithmic interaction. This is not at all the case for A = 1 and 10 (see FigsjHl^a) and (b)), 
but for A = 50 it does. Therefore, transition to the logarithmic-type of interaction takes 
place between A = 10 and 50. 

It is also interesting to look at the eigenmodes corresponding to the highest frequency for 

= 20 (see Figs.E^a-c)). Now the modes for A = 1, 10 and 50 look all similar. However, they 
differ from the breathing mode: the middle and outer ring move in the opposite direction. 
From the spectrum one can observe that for > 10 there are two highest frequencies which 
lie very close together. These highest but one frequencies are therefore plot in Fig. IHfd-f). 
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Notice that for A = 10 and 50 this mode indeed corresponds with the breathing mode. 

VI. CONCLUSIONS 

We investigated the ground state structure and the spectral properties of a classical 
2D cluster with a finite number of particles interacting through a repulsive logarithmic 
potential. The logarithmic potential is the Coulomb potential in a 2D world. It was shown 
that for a hard wall confinement all particles are situated on one ring at the boundary of the 
potential, as a consequence of Earnshaw's theorem in 2D. Multi-ring configurations exist for 
soft confinement potentials in which case only the external confinement delivers the restoring 
force. The particles are uniformly distributed for a parabolic confinement. 

We also found that, independently of the number of particles in the cluster, all eigen- 
frequencies lie below the frequency of a particular mode, namely the breathing-like mode. 
This is again a consequence of Earnshaw's theorem because adding extra particles does not 
create extra minima in the potential nor makes it more complex. In the case of a parabolic 
confinement the frequency of this breathing mode is independent from N, which results in 
an upper bound for the frequencies. 

It is also shown that for vortices in a type II superconductor (for which the penetration 
length A ^ coherence length) the transition to the logarithmic-type of interaction takes 
place between 10 < A < 50. It is not always the highest mode anymore which corresponds 
with the breathing mode, it can also be the highest but one mode. 
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